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,!_H | 1 INTRODUCTION 



We present low resolution (R ~ 450) -fCband spectroscopy for 16 of the 43 cir- 
cumnuclear star-forming knots in M100 identified by Ryder & Knapen (1999). We 
compare our measurements of equivalent widths for the Br7 emission line and CO 
2.29 /im absorption band in each knot with the predictions of starburst models from 
the literature, and derive ages and burst parameters for the knots. The majority of 
these knots are best explained by the result of short, localised bursts of star formation 
between 8 and 10 Myr ago. By examining both radial and azimuthal trends in the 
age distribution, we present a case for sequential triggering of star formation, most 
likely due to the action of a large-scale shock. In an appendix, we draw attention to 
the fact that the growth in the CO spectroscopic index with decreasing temperature 
in supergiant stars is not as regular as is commonly assumed. 

Key words: galaxies: evolution - galaxies: individual (M100, NGC 4321) - galaxies: 
kinematics and dynamics - galaxies: structure - infrared: galaxies. 
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Nowhere is the interplay between galaxy dynamics and star 
formation more important than in the circumnuclear regions 
(CNRs) of barred spiral galaxies. Where star formation is 
observed to occur, it is frequently organized into a tightly- 
wound spiral, or even a complete ring, on the scale of a kilo- 
parsec or more. Based on numerical simulations of gas flows 
in a barred potential, it is thought that these features are 
associated with dynamical resonances, specifically the Inner 
Lindblad Resonances (ILRs). The resulting concentration of 
gas in a narrow region leads eventually to the triggering of 
star formation in a host of individual 'hot-spots' or 'knots'. 
However much of the detailed specifics, such as whether star 
formation commences in all such knots at the same time, and 
in the same manner, remain a mystery. 

In an effort to address such issues, we recently obtained 
(Ryder & Knapen 1999, hereafter Paper I) high-resolution 
images in the J, H, and K bands with IRCAM3 on the 
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United Kingdom Infrared Telescope (UKIRT) of the circum- 
nuclear region of M100 (NGC 4321), the brightest spiral 
galaxy in the Virgo cluster. In conditions of < 0.4 arcsec 
seeing, we were able to identify 43 compact knots in the 
K image, and to determine magnitudes and colours for 41 
of these. Spectroscopic observations with CGS4 on UKIRT, 
with the slit position angle set to cross the nucleus and 
the two large-scale hot-spots Kl and K2 identified previ- 
ously by Knapen et al. (1995a), revealed the presence of 
weak Br7 and H2 1-0 S(l) emission lines, as well as deep 
CO (2-0) absorption bands in both Kl and K2. A compari- 
son with models for starburst galaxies by Puxley, Doyon & 
Ward (1997; hereafter PDW) indicated that these knots were 
formed between 15 and 25 Myr ago in a burst of star forma- 
tion which decayed exponentially with a 5 Myr timescale. 
However, with spectroscopic data for only 2 such regions, it 
is difficult to draw useful conclusions about the nature of the 
circumnuclear starburst in M100. Since multi-object, near- 
infrared spectrographs are still under development, long-slit 
spectroscopic observations like those presented here are of 
necessity somewhat selective, but we have been able to se- 
cure spectra for more than one-third of the compact knots 
discovered in Paper I. 

It is possible to draw some conclusions regarding the 
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circumnuclear star formation history of spiral galaxies from 
imaging in broad-band filters. However, broad-band colours 
are extremely sensitive to the amount of reddening assumed, 
as well as to any emission from hot dust. Narrow-band 
imaging of emission lines, including the use of Fabry-Perot 
ctalons (e.g., Reunanen et al. 2000; Kotilainen et al. 2000), 
provides a potentially more sensitive age diagnostic. Models 
have demonstrated that the strength of the Br7 emission 
line equivalent width begins to decline ~ 3 Myr after the 
onset of a burst, but the rate of decline will depend on the 
decay rate of the starburst, the IMF, the metallicity, etc. 

Over the same period, the equivalent width of the CO 
(2-0) bands will grow to a maximum, and then decline, lead- 
ing to an ambiguity in age which is also model-dependent. 
Recently, the very usefulness of the CO band strengths as 
an age- indicator for star formation has come under renewed 
scrutiny (Origlia et al. 1999; Origlia & Oliva 2000), princi- 
pally because of still-incomplete handling of the AGB phase 
of red supergiant evolution by the evolutionary models, as 
well as the competing effects of stellar effective tempera- 
ture, surface gravity, metallicity, and microturbulent veloc- 
ity. Nevertheless, for starburst ages greater than 5 and up 
to 100 Myr, we aim to demonstrate that the combination of 
Br7 and CO (2-0) equivalent width measurements can place 
quite useful constraints on the available models, and serve 
as a powerful relative age discriminant. 

In this paper, we present new low-resolution K-hanA 
spectra for many more of the circumnuclear star-forming 
knots in M100, and determine ages and burst parameters 
for them. We begin by describing our spectroscopic obser- 
vations and measurements in Sections ^ and |j| In Section 
we compare our observed line equivalent widths with star- 
burst models, and use the derived ages to search for radial 
and azimuthal trends. Our conclusions are contained in Sec- 
tion ^, and in an Appendix, we highlight some aspects of 
the CO spectroscopic index in nearby stars that are still not 
generally appreciated. 



2 OBSERVATIONS AND DATA REDUCTION 

Longslit spectroscopy of the central region of M100 was car- 
ried out over three nights from 1999 April 1-3 UT. Only 
limited data was acquired the first night due to heavy cir- 
rus; conditions were photometric throughout the remaining 
two nights. The common-user near-IR spectrograph CGS4 
(Mountain et al. 1990) on UKIRT with its 0.61 arcsec pixel 
scale, a 1.2 arcsec slit, and the 40 line mm -1 grating in first 
order was employed, which delivers complete spectral cov- 
erage from 1.85-2.45 /xm at a resolving power of 450. Using 
the images presented in Paper I, a total of 4 slit position 
angles and nuclear offsets were selected, covering at least 

3 knots each; these are marked on Figure 

After taking the observations presented in Paper I, it 
became apparent that the infrared-emitting central region 
of Ml 00 fills no more than one-third of the usable area of 
the detector array, and therefore, nodding to blank sky was 
not really necessary. For the observations presented here, sky 
subtraction was performed by sliding the objects of interest 
a total of 60 rows (36.6 arcsec) along the slit in between 
'object' and 'sky' exposures. After setting the slit to the 
desired position angle, a manual search for the K-h&nd nu- 



clear peak was conducted, following which accurate offsets 
were applied by use of the UKIRT crosshead to place the 
slit right at the locations marked on Fig. [lj All observations 
of M100 were bracketed by similar sequences on the nearby 
A4 V star BS 4632, to enable removal of some of the telluric 
absorption features. 

Data reduction was accomplished using a combination 
of tasks within the Starlink CGS4dr and NOAO iraf^| pack- 
ages. After flatfielding, the difference of each object-sky pair 
was co-added into groups, with each grouped image repre- 
senting ~ 3000 sec of on-source integration in the case of 
the M100 spectra, or ~ 300 sec for BS 4632. Any residual 
sky absorption/emission was removed from this image by 
fitting a low-order polynomial to galaxy-free regions within 
each row of the spatial axis. The high signal-to-noise of the 
stellar continuum in BS 4632 was used to define the width 
and curvature of suitable extraction apertures for both the 
'positive' and 'negative' beams. After extracting matching 
apertures from an argon arc image, the spectra from both 
beams were wavelength calibrated separately to an accuracy 
of 0.2 nm. The spectrum of the negative beam was inverted, 
and then co-added in wavelength space to that of the posi- 
tive beam. The strong Br7 absorption intrinsic to BS 4632 
was interpolated over, but this could not be done for Paa 
or for Br5 due to their location near the trough of an H2O 
absorption band. 

Processing of the M100 knot spectra proceeded in a 
similar fashion, except that multiple apertures needed to be 
defined for each slit position. This was done in consulta- 
tion with Fig. |lj but the exact aperture centres and widths 
were defined independently for each group. After extrac- 
tion, wavelength calibration, and co-addition, the spectrum 
of each knot was divided by the BS 4632 spectrum which 
best matched the airmass of observation, and then multi- 
plied by a blackbody curve corresponding to the temper- 
ature (8380 K) and flux density (K = 6.13) of BS 4632. 
Table hi summarises the observational parameters for each 
slit setting, including all of the knots that each slit was ex- 
pected to intersect to some degree. The numbering sequence 
for the knots is as given in Table 1 of Paper I. 

In order to define a consistent continuum right out into 
the region of extensive CO absorption beyond 2.3 /im, we fit- 
ted a power-law of the form F\ oc X 13 to featureless sections 
of the spectrum near 2.11 and 2.27 /im (rest wavelength). 
After normalising each spectrum by this fit, the equivalent 
widths of the redshifted Paa, Br7, and H2 lines were mea- 
sured using the SPLOT task in iraf, together with two de- 
terminations of the spectroscopic index CO sp . The original 
definition of a CO spectroscopic index by Doyon, Joseph & 
Wright (1994a; hereafter DJW) has the form 

CO? p JW = -2.51og(i? 2 . 36 > (1) 
where -R2.36 is the mean normalised intensity between 2.31 



t IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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and 2.40 /im in the rest frame of the galaxy (in the case 
of M100, this wavelength interval corresponds to 2.322 — 
2.412 /im). Although the output of both the DJW and the 
Leitherer et al. (1999; hereafter SB99) models use this def- 
inition of CO sp , PDW advocated the adoption of a nar- 
rower wavelength interval over which to integrate, namely 
2.2931 — 2.3200 /jm, which covers just the region from the 
band head to near the end of the 2-0 band. After a com- 
prehensive analysis, this range was found to offer the best 
sensitivity to changes in luminosity class, and has also been 
endorsed by Hill et al. (1999). From equation 7 of PDW, we 
get the following relation between their definition of CO sp 
and the equivalent width of CO (nm), measured over the 
rest wavelength interval 2.2931 — 2.3200 ^m: 



COf 



-2.5 log 



16.8 - EW(CO,nm) 



(2) 



coupled with the need to then apply an empirical transfor- 
mation, leads us to favour the use of CO^, JW for comparison 
with the models. While the PDW definition undoubtedly is 
more robust than even narrower wavelength regions, and is 
more sensitive to changes in the stellar population, there ap- 
pears to be little advantage to using it at our resolution. Our 
conclusions in the next section regarding the relative ages of 
the star-forming knots in M100 are unaffected by the choice 
of CO spectroscopic index. 

in the range 0.20-0.28, 

we infer from Table 



Al 



that the stellar populations of these 
knots are currently dominated by giants and/or supergiants 
of type K5 or later. We do not observe values of CO^, JW 
quite as high as was measured for the 'hotspots' Kl and K2 
in Paper I, which we attribute to the fact that the original 
spectroscopy did not target any particular knot, and that 
the continuum fitting was less rigorous. 



At the redshift of M100 (1571 km s" 1 ; Knapen et al. 1993), 
this wavelength interval becomes 2.305 — 2.332 /im. Note 
that PDW give r\ = 16.6, but our analysis in the Appendix 
shows that r\ = 16.9 is more appropriate at our resolution. 



3 RESULTS 

Figure || shows the spectra of the combined data for three 
of the knots, on which the main spectral features that have 
been measured are identified. Even at this resolution, the 
redshift of M100 is just enough to separate the Paa emis- 
sion line from the anomalous bump at 1.876 /im, which is a 
residual of the Paa absorption in BS 4632. The (almost flat) 
continuum slope found in most of the knot spectra helps 
greatly in minimising the error that comes from fitting a 
power-law to the continuum and extrapolating it across the 
CO band. 

Table ^| is a compilation of the measured equivalent 
widths and continuum power-law slopes (3 for each of the 
knots which could be clearly distinguished in the spec- 
tra. Because of the tendency for COfjP w to underestimate 
CO^, JW as noted in the Appendix, the former has been cal- 
culated using r\ = 16.9 in equation |2| In order to gain a real- 
istic estimate of the uncertainties, the equivalent width mea- 
surements were performed on each of the 3 or 4 individually- 
reduced groups (with exposure times of 2400-3800 sec), as 
well as the co-added result. The mean equivalent width 
adopted is that measured from the co-added spectra, while 
the error bars are based on the dispersion in the individ- 
ual measurements. As expected, this shows that the uncer- 
tainties are dominated as much by the noise statistics and 
the subjective choice of extraction aperture as by the mea- 
surement errors and continuum fitting. Note that we have 
two independent observations and measurements of knot 28 
(with slit 1 and, less optimally, slit 2), and except for the 
equivalent width of H2 , the results are consistent within the 
error bars. 

It is also apparent that the error bars on measurements 
over the smaller PDW wavelength interval (their "extended" 
region) are significantly larger than the error bars on mea- 
surements over the full DJW range. These larger error bars, 



4 DISCUSSION 

4.1 The Ages of the Starburst Knots 

Figure |^ presents our results as a plot of CO spectroscopic 
index vs. the Br7 equivalent width for the 16 knots (the 
two independent observations of knot 28 have been aver- 
aged), together with the loci of evolutionary synthesis mod- 
els for starbursts as calculated by DJW, and used by PDW 
in age-dating the nuclear starburst in M83. There appears 
to be a well-defined sequence in the line strengths among 
the bulk of the knots, which parallels quite closely the star- 
burst model with the shortest burst timescale r (i.e. the 
"quasi-instantaneous" burst). Reading off the fiducial marks 
on this model indicates ages for these knots of between 8 and 
10 Myr. 

There is of course a selection effect here: the Br7 equiva- 
lent width must be at least 0.2 nm to be detectable, so we can 
only put a lower age limit on knot 13. Similarly, there could 
be circumnuclear star-forming knots younger than 5 Myr 
around M100 which would be missed by our NIR imaging 
and spectroscopy, since there has not been sufficient time for 
the most massive stars to evolve to red supergiants. Never- 
theless the observed upper envelope of Br7 equivalent widths 
just below 1 nm suggests that there was an extended pause 
in star formation events ~ 8 Myr (on this model) ago. Fur- 
thermore, any bursts decaying at a rate slower than 5 Myr, 
and aged between 10 and 30 Myr ought also to be detectable, 
but none are seen. 

We have also plotted our results on the closest equiv- 
alent SB99 models (Figure ^J). The agreement between ob- 
servations and the models is not as good as for the DJW 
models, but the two evolutionary tracks shown are the only 
two that come anywhere close to achieving the required val- 
ues of CO™. The first model is for an instantaneous burst 
of star formation (ISF), in which the IMF has a Salpeter-like 
form (slope a = 2.35 and upper mass cutoff of 100 M©) and 
a metallicity of 2x solar (Z — 0.04). The alternative model 
is one which has a constant star formation rate (CSF) of 
1 Mq yr _1 , a steeper IMF (a = 3.30 and upper mass cutoff 
of 100 Mq) and a metallicity of twice-solar. 

A high central metallicity in M100 is supported by 
H 11 region spectroscopic abundance measurements by Skill- 
man et al. (1996), who found 12 + log(0/H) ~ 9.3 near 
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the nucleus, or over 2.5 x solar metallicity. The tendency of 
existing models to predict CO spectroscopic indices which 
are too weak is a known problem (e.g., Devost & Origlia 
1998; Origlia et al. 1999), and appears to be related to un- 
certainties in the effective temperatures and microturbulent 
velocities of the red supergiants, as well as the fraction of 
time these stars spend as blue supergiants during their core 
helium burning phase. Nevertheless, the evolutionary trend 
of the ISF model parallels quite closely the distribution of 
the knots, and even the age range of 6.5 — 7.0 Myr is not 
inconsistent with that indicated by the DJW models, con- 
sidering that the ISF model has r < 0.1 Myr. 

Although the (quasi-)instantaneous burst models fit the 
observed trend of knot line indices quite well, both Fig- 
ures ^ and § leave open the possibility that a few of the 
knots could equally well just be due to steady (or only 
slowly declining) continuous star formation for the past 
50-100 Myr. The same CSF model plotted in Fig. § pre- 
dicts that in this age range, the knots should have colours 
(J-H) = 0.48 ± 0.02 and (H - K) = 0.24 ± 0.02, while the 
ISF model in the same figure predicts (J — H) — 0.46 ± 0.05 
and [H - K) = 0.22 ± 0.04 for the age range 6.6-7.0 Myr. 
Each of these cases is fully consistent with the median of 
the dereddened knot colours, as plotted in Fig. 2 of Paper I, 
emphasising how ineffective broadband near-infrared colours 
are in distinguishing bu rsts fr om continuous star formation. 
As we show in Section 4.2.3, even knowledge of the H2 1— 



S(l) to Bi'7 ratio still does not allow for an unambiguous 
distinction between instantaneous burst and continuous star 
formation models. 



4.2 Other Hydrogen Lines 

4.2.1 Pact 

As Fig. ^ shows, the Paa line at the redshift of M100 
(1571 km s _1 ) is clearly resolved from the rest-frame Paa 
absorption in BS 4632. Thanks to the low water vapour col- 
umn over Mauna Kea, measurements of Paa are quite fea- 
sible despite the low atmospheric transmission at this wave- 
length, although the continuum level in this region is not 
as well fit by the power law. For this reason, the equivalent 
widths in Table [2 have been measured with respect to the 
base of the line profile, rather than the continuum as defined 
by an extrapolation of the power law with index /3. 

In Figure ^| we have plotted the equivalent width of the 
Paa line against that of the Br7 line. For Case B recombi- 
nation of hydrogen at T = 10 4 K (e.g. Table 4.2 of Oster- 
brock 1989), we would expect a ratio J(Pa)/7(Br7) = 12.52, 
as indicated by the dashed line (more metal-rich gas will 
be cooler, but even at T = 5 x 10 3 K, the ratio is only 
slightly higher at 12.86). The fact that the Paa line is con- 
sistently weaker than Br7 would predict is usually attributed 
to reddening. However, the mean amount of extinction re- 
quired to explain this difference (adopting the interstellar 
extinction curve from Cardelli, Clayton & Mathis 1989) is 
Ay ~ 16 mag, or Ak ~ 1.8 mag. This is a factor of 4 larger 
than the maximum extinction that was derived from the 
near-IR colour excesses of the knots in Paper I. 

The equivalent width of an emission line ought to be 
independent of reddening, but this assumes that both the 
continuum flux and the line emission at a given wavelength 



originate from the same physical region; if the red super- 
giants had cleared much of the dust that shrouded their 
birth, then the stellar and nebular fluxes could experience 
quite different amounts of reddening. Alternatively, our cor- 
rections for Paa absorption in the standard star, and for 
H2O absorption in our atmosphere, may be inadequate. An 
added complication is the presence of extended, diffuse Paa 
emission in the nuclei of many early-type galaxies (Boker et 
al. 1999; Alonso-Herrero et al. 2001). Clearly, despite being 
an order of magnitude stronger than Br7, there are still a 
number of difficulties with using ground-based observations 
of the Paa line as a star formation and reddening diagnostic. 



1 2. 2 Ha 

The Ha narrow-band imaging and Fabry-Perot kinematics 
presented by Knapen et al. (1995a; 2000) reveal four main 
complexes of H 11 regions forming an almost complete nuclear 
ring, coincident with that outlined by the Ji'-band knots. 
Two of these complexes are associated with the major 'hot- 
spots' Kl and K2 at either end of the bar. On the basis 
of reduced dust obscuration, they suggested that Kl ought 
to be slightly older than K2, and indeed, our age analy- 
sis confirms that the aggregate of knots which make up Kl 
(knots 28-31) is slightly older than that of K2 (knots 19 and 
20). Complexes Ha3 and Ha4, along the bar minor axis, are 
brighter than Kl and K2 in Ha, but dimmer in K, consistent 
with their being less evolved. As discussed in Section ^l], the 
youngest knots are found near Ha3 and Ha4, providing fur- 
ther support for the qualitative age sequence of Knapen et 
al. 



4-2.3 Molecular hydrogen 

Molecular hydrogen emission in the H2 1-0 S(l) line at 
a rest wavelength of 2.122 /jm is commonly observed in 
star-forming environments, but the nature of its excita- 
tion (shocks vs. UV fluorescence) is still somewhat ambigu- 
ous (e.g., Doyon, Wright & Joseph 1994b; Puxley, Ram- 
say Howat & Mountain 2000). When we plot the H2 1- 
S(l) line strengths against Br7 (Figure m), the two are 
clearly correlated (just as was demonstrated for the nuclei of 
a large sample of star-forming galaxies by Puxley, Hawarden 
& Mountain 1990), suggesting a close link between molecu- 
lar hydrogen line emission and star formation in these knots. 
Table |i| lists the ratios of the H2 1-0 S(l) to Br7 equivalent 
widths, and these are predominantly in the range 0.4 — 0.9, 
very similar to the Puxley et al. (1990) ratios. 

Using the same modeling framework as used for the 
evolutionary tracks in Fig. ^, Doyon et al. (1994b) showed 
that the S(l)/Br7 ratio would not exceed 0.4 until at least 
60 Myr after the onset of the burst, for r = 20 Myr. The 
short (t — 1 Myr) bursts favoured by Fig. ^ might be ex- 
pected to achieve the observed S(l)/Br7 ratios in Table ^ 
even earlier than 60 Myr. Their model assumes however that 
outflows from young stellar objects (YSOs) and from su- 
pernova remnants (SNRs) produce all of the H2 line flux, 
whereas both UV fluorescence of photodissociation regions 
by massive stars, and large-scale shocks from the action of 
the spiral density wave, may also be significant sources in 
the centre of M100. 
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For comparison, their model with continuous star for- 
mation (like the CSF model of SB99 in Fig. |) would barely 
exceed a ratio of 0.15, even after 100 Myr of evolution. Vanzi 
& Rieke (1997) measured ratios of this order in a sample 
of blue dwarf galaxies which, while representing perhaps 
the purest form of a starburst, may also be subject to the 
unknown effects of low metallicity. Thus, even a combined 
knowledge of the knot near-IR colours, their Br7, H2 and 
CO equivalent widths is still not sufficient to definitively ex- 
clude the possibility that at least some of the circumnuclear 
knots in M100 have been continuously forming stars for as 
long as 100 Myr. 

4.3 Age trends 

Table ^| lists the ages derived for each of the knots which lie 
on, or very near to the r = 1 Myr locus of DJW in Fig |S| 
Four of the knots (2, 13, 28, and 38) lie well off this se- 
quence, and it would appear that the bursts powering them 
are declining at a different rate from all the rest. Clearly, 
ages derived from the SB99 ISF model in Fig. ^ would be 
slightly younger, but the relative age sequence will still be 
the same. The location of each knot within the deprojected 
disk of M100 has been calculated using the disk orienta- 
tion parameters [i = 27°, (j> = 153°) derived by Knapen 
et al. (1993) from their Hi kinematic analysis. In Figure 
we plot these ages against both radial distance and posi- 
tion angle in the plane of M100. In this sense, position angle 
increases counter-clockwise, with the origin lying midway 
between knots 20 and 37 in Fig |jj 

We consider first the azimuthal age distribution. The 
knots at either end of the bar (19, 20 and 29) are all about 
the same age, ~ 8.8 Myr. The oldest knots (30, 31, 34, and 
35) are located just south of the western end of the inner 
bar, while the youngest knots (9, 37 and 39) lie close to the 
bar minor axis, leading (or lagging) the oldest knots by 90° . 
Indeed, beginning at position angle 70° , there is the sugges- 
tion of an age sequence spanning at least a Myr as one circles 
the nuclear 'ring'. A simple linear least squares fit (Press et 
al. 1992) to the age distribution, starting at position an- 
gle 70°, yields a Pearson's r value of —0.75 and Student's t 
probability of 0.005, the small value of which implies a sig- 
nificant correlation between position angle around the ring, 
and burst age. However, one could just as easily point to the 
appearance of two age minima (at position angles 50° and 
230°), but our limited sampling prevents us from identify- 
ing if there is a corresponding second age maximum near 
300°. Nevertheless, Fig. Q provides the first evidence that 
rather than being a stochastic process, circumnuclear star 
formation may be sequentially triggered. 

At the radial distance of these knots, the orbital pe- 
riod given by the rotation curve of Knapen et al. (1993) is 
~ 20 Myr, an order of magnitude larger than the observed 
age spread. Similarly, the bar pattern speed (determined in- 
dependently by Knapen et al. (1995b) and by Wada et al. 
(1998)), Q p ~ 65 km s" 1 kpc -1 leads to an even longer 
timescale. The maximum streaming velocities observed in 
Ha and CO by Knapen et al. (2000) around the circum- 
nuclear region are of order 40 km s _1 , too slow to causally 
connect the oldest and youngest star-forming knots. Thus, 
it cannot be the rotation of the bar, or bulk motion of the 
gas, which directly triggers the star formation. 



The detailed nuclear morphology of M100, as pieced 
together from optical, near-infrared, Ha, and CO imaging, 
is analysed in great detail by Knapen et al. (1995a, b; 2000). 
They constructed a numerical model including gas, stars, 
and star formation, which successfully replicated many of 
the main features, including (1) the 1 arcmin stellar bar; 
(2) a pair of offset curved shocks, and streaming motions 
along the nuclear 'ring', both signatures of a global density 
wave induced by this bar; and (3) four distinct compression 
zones, two of which are close to the bar minor axis (where 
local maxima in the CO and Ha emission are observed), 
and two near the ends of the bar where the spiral armlets 
are observed to switch from leading to trailing. As shown in 
fig. 15 of Knapen et al. (1995b), star formation at the ends 
of the bar in their model is seen to precede that along the 
minor axis (just as Fig. |?] suggests) , but again the timescales 
involved are rather longer than the observed age spread. 

One scenario worth considering would be an outward 
propagating shock-wave, which reaches the ends of the bar 
first, then sets off the knots along the bar minor axis a short 
while later, resulting in a radial age gradient. As Fig. ^ 
shows, the oldest knots do tend to be at larger radii, but 
then so too is the youngest, so the evidence is less com- 
pelling (for a linear least squares fit, r — 0.22 and t = 0.5). 
What makes this model attractive though, is that a shock 
wave velocity of only ~ 200 km s _1 is required to explain the 
observed age spread, compared with up to 1000 km s _1 for 
a density wave-driven shock propagating around the ring. 
In addition, a one-off outburst such as this makes it easier 
to account for the apparent recent pause in star formation 
events, whereas a persistent phenomenon (such as a spiral 
density wave for instance) would be expected to have pro- 
duced a much larger spread in burst ages than is observed 
in M100. 



5 CONCLUSIONS 

New longslit i("-band spectroscopy of the circumnuclear 
region of M100 has been obtained, targeting specifically 
the compact 'knots' identified from high-resolution near-IR 
imaging in Paper I. We have shown that even at compara- 
tively low resolution (R ~ 450), it is possible to discern age 
differences of as little as 0.2 Myr in circumnuclear star form- 
ing regions by comparing the equivalent widths of the Br7 
emission line with the CO 2.29 /im absorption band. While 
the primary age sequence diagnostic is the CO spectroscopic 
index, complementary measurements of the Br7 and H2 1— 
S(l) lines allows us to distinguish short-duration, rapidly 
declining bursts from continuous star formation. Compari- 
son with starburst models from the literature indicates that 
the bulk of the circumnuclear star formation events in M100 
are best accounted for by the former, with decay timescales 
of ~ 1 Myr. 

We find the strongest evidence yet that circumnuclear 
star formation in M100 is sequentially triggered in an az- 
imuthal or radial sense (or perhaps both). While the ab- 
solute ages of the knots are highly model-dependent, the 
age spread is unambiguous. The youngest and oldest knots 
are quite spatially distinct, but the physical separations and 
timescales involved are such that only shocks with veloci- 
ties of a few hundreds of kilometres per second can causally 
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connect them. The presence of curved dust lanes bisecting 
the circumnuclear ring verifies the existence of large-scale 
shocks in this region, most likely driven by a bar-induced 
spiral density wave. However, the symmetry of the age dis- 
tribution is such that a simple, radially-propagating shock 
may just as easily explain the situation, though no source 
can yet be identified. 

To date, we have secured spectra for 16 of the 43 knots 
identified in Paper I, and of these, only 12 are suitable for 
constraining the age distribution. To confirm our suspicion 
of an azimuthal age gradient (or even a cycle) in M100 will 
require better sampling of some of the fainter knots in M100 
(something for which the new generation of near-infrared 
multi-object spectrographs and integral-field units will be 
ideally suited). Few galaxies have had their circumnuclear 
kinematics and star formation studied in as much detail as 
M100 has, but more are clearly needed if we are to under- 
stand the role of the bar in influencing both. 
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Table 1. CGS4 Slit Parameters for M100 



Slit 


P.A. 


Knots 


Total Exp. (sec) 


1 


169.6 


28, 29, 34, 35, 38 


10080 


2 


17.9 


2, 28, 30, 31 


14880" 


3 


128.3 


9, 10, 11, 13, 14 


7200 


4 


60.9 


18. 19, 20, 37, 39 


8640 



"Includes 5520 sec through thick cirrus. 



Table 3. Knot Locations and r = 1 Myr Burst Ages 



Knot 


R 


e 


Age 




(arcsec) 


(°) 


(Myr) 


2 


13.0 


190 




9 


6.9 


232 


8.4 ±0.1 


10 


6.6 


249 


8.8 ±0.2 


11 


7.0 


263 


8.8 ±0.2 


13 


8.6 


286 




19 


7.9 


319 


9.0 ±0.1 


20 


7.2 


326 


8.8 ±0.2 


28 


8.2 


149 




29 


7.6 


142 


8.8 ±0.2 


30 


8.7 


131 


9.6 ±0.4 


31 


9.3 


118 


9.6 ±0.2 


34 


5.9 


116 


9.2 ±0.2 


35 


7.0 


82 


9.3 ±0.2 


37 


6.5 


27 


8.6 ±0.1 


38 


8.6 


59 




39 


9.3 


47 


8.3 ±0.2 
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Knot 





-EW(Pa) 
(nm) 


-EW(H 2 1-0 S(l)) -EW(Br 7 ) 
(nm) ( nm ) 


S(l)/Br 7 


CO? p ' w 
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Slit 3 
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10 
11 
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0.46 ±0.12 
< 0.21 


0.53 ±0.11 

< 0.07 
1.46 ± 0.51 

< 1.1 


0.20 ±0.01 
0.23 ±0.01 
0.22 ±0.01 
0.26 ±0.01 
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Slit 1 



Figure 1. Location of the CGS4 slit positions used, together with identifications of the knots observed, marked on the _R"-band image 
of the nuclear region of MfOO from Ryder & Knapen (1999). The 1.2 arcsec slit width is shown to scale. 
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Figure 2. Reduced spectra for three of the knots in M100, illustrating the relative changes in line strengths going from a relatively young 
knot (Knot 9, age ~ 8.4 Myr) to one of the oldest knots (Knot 30, age ~ 9.6 Myr), via one of intermediate age (Knot 20, age ~ 8.8 Myr). 
The main spectral features are marked, as well as the two wavelength regions over which the CO equivalent width was measured: {upper 
line) the range defined by DJW, and (lower line) the narrower range adopted by PDW. The bump on the blue wing of the Pact lines is 
a residual of the Pact 1.876 £im absorption intrinsic to the telluric standard star. The residual due to Br<5 at 1.945 (im has been patched 
over in these spectra. 
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Figure 3. Observed CO spectroscopic indices plotted against the Br7 equivalent widths for 16 of the star-forming knots in M100. The 
two independent observations of knot 28 in Table ^ have been combined with a weighted average, to give the point at (—0.73, 0.23). Also 
shown in this plot are the evolutionary tracks for three starburst models as presented in Fig. 3 of PDW, with exponentially-decreasing 
star formation timescales r = 1, 5, and 20 X 10 6 yr. The points marked on the tracks indicate the time taken in Myr to reach that point 
from the onset of the burst. 
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log Bry EW (nm) 



Figure 4. As for Figure H, but this time using two evolutionary tracks from Leitherer et al. (1999) models. The track marked "ISF" 
is for a Z = 0.04, a = 2.35, M up = 100 Mq instantaneous burst model, while the one marked "CSF" is for a Z = 0.04, a = 3.30, 
M up = 100 Mq model with continuous star formation at a rate of 1 Mq yr — 1 . The points marked on the tracks indicate the time taken 
in Myr to reach that point from the onset of star formation. 
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Figure 5. The equivalent width of the Paa emission line from each of the knots plotted against the Br7 equivalent width. The relationshi 
expected from Case B recombination line theory is shown by the dashed line. The two independent observations of knot 28 in Table 
have been combined with a weighted average, to give the point at (—0.73, 0.29). 
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Figure 6. The equivalent width of the H2 1-0 S(l) emission line from each of the knots plotted against the Br7 equivalent width. The 
two independent observations of knot 28 in Table have been combined with a weighted average, to give the point at (—0.73, —0.48). 
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Figure 7. The derived ages for the 12 knots in M100 which appear to best follow the evolutionary trend of the r = 1 Myr burst locus 
in Fig. |i| plotted against deprojected radius (top) and azimuth (bottom) in the disk of M100. The zero-point of the position angle scale 
corresponds to the kinematical major axis <f> = 153° (roughly halfway between knots 20 and 37 in Fig. lh, and position angle increases 
going counter-clockwise. 
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APPENDIX A: THE CO INDEX IN STARS 

Since many of the studies to date of the strength of CO 
absorption in starburst systems have been with resolving 
powers at least twice that which we have used here, we felt 
it was important to verify that we would not lose any sen- 
sitivity to stellar effective temperature or luminosity class. 
To this end, we set out to observe a variety of late-type 
stars covering nearly the full range of spectral type later 
than F5, with the same observational setup as for the obser- 
vations of M100. Observing time constraints, coupled with 
the difficulty in finding cool stars with well-defined spec- 
tral types which were not so bright as to saturate the de- 
tector, precluded us from going below T e g = 3400 K. Ta- 
ble Al gives the observed CO sp values for both definitions 
(equations |l| and ^), as well as the spectral classification 
from Hofneit (1982), the derived stellar effective tempera- 
ture from Schmidt-Kaler (1982), and the continuum power- 
law slope /3. Figure Al presents the results for our stellar 
dataset, along with the relations for dwarfs, giants, and su- 
pergiants empirically-determined by DJW from the stellar 
atlas of Kleinmann & Hall (1986). 

Although our measurements match the expected trends 
quite well for the hotter stars, there is considerably more 
scatter than for the equivalent plot in DJW (their Fig. 7). 
In order to assess how much of this scatter is due to our res- 
olution, and how much may be intrinsic to the stars them- 
selves, we have made use of the if-band stellar spectral li- 
brary of Wallace & Hinkle (1997). We measured CO° JW 
in all of their luminosity class I, III, and V stars with 
3000 < T c ff < 6000 K that w ere n ot also in the Kleinmann & 



Table Al. Stellar Parameters and CO Spectroscopic Indices 



Hall (1986) study. As Figure A2 shows, even at higher reso- 
lution (R ~ 3000) there is indeed somewhat more scatter in 
CO spectroscopic index with type and luminosity class than 
one might conclude from the analysis of DJW, particularly 
amongst the bright supergiants (a similarly enhanced degree 
of scatter amongst the supergiants is exhibited by the inde- 
pendent sample of Forster Schreiber 2000) . While this does 
provide some reassurance that lower resolution does not sig- 
nificantly compromise our ability to measure real changes in 
CO^,' TW with changing stellar population, the large scatter 
(up to 0.1 dex) about the mean relation defined for super- 
giants by DJW ought to be of concern. The supergiants, 
when present, will dominate CO^, JW , and the predictions of 
the models assume a much tighter relation between CO^, JW 
and T e fr like the one found by DJW. 

One aspect that may be resolution-dependent is the 
conversion from CO equivalent width over the narrow wave- 
length interval used by PDW, and the equivalent CO spec- 
troscopic index (equation [2J). Lowering the resolution will 
tend to broaden the CO bandhead, and the wavelength lim- 
its may no longer encompass the same extent of the line pro- 
file. The much larger wavelength interval used by DJW will 
be less susceptible to this effect. Indeed, when we plot CO sp 
using both definitions for our own ob serv ations of field stars, 
as well as all the M100 knots (Figure A3), it is apparent that 



CO^ p J w is consistently less than COfp ^ by ~ 0.02 dex when 
equation ^ is used with 77 = 16.6. Hill et al. (1999) found a 
similar tendency to underestimate CO^ p DW in their R ~ 400 
spectra. We can compensate for this by using 77 = 16.9 in 
equation |^ to restore good agreement between the two in- 
dices. 



Star 


Type 


P 


Toff 


CO? p JW 


COP p DW 


BS 7503 


G1.5 Vb 


-1.5 


5900 


0.00 


-0.01 


BS 6697 


G2 V 


-1.4 


5860 


0.01 


-0.01 


BS 6538 


G5 V 


-1.4 


5770 


0.01 


-0.01 


BS 6301 


K0 V 


-1.3 


5250 


0.04 


0.02 


BS 6806 


K2 V 


-1.6 


4900 


0.03 


0.00 


BS 8085 


K5 V 


-1.2 


4350 


0.06 


0.05 


BS 8086 


K7 V 


-1.1 


4060 


0.06 


0.06 


BS 6531 


F6 III 


-1.7 


6360 


0.00 


-0.02 


BS 6466 


GO III 


-1.6 


5850 


0.04 


0.02 


BS 6239 


G5 III 


-1.5 


5150 


0.04 


0.01 


BS 6287 


G8 III 


-1.9 


4900 


0.04 


0.02 


BS 6307 


K0 III 


-1.4 


4750 


0.07 


0.04 


BS 6358 


K5 III 


-2.0 


3950 


0.14 


0.12 


BS 6159 


K7 III 


-2.3 


3850 


0.17 


0.15 


BS 7244 


M0 III 


-2.4 


3800 


0.18 


0.12 


BS 6834 


M4 Illab 


-2.1 


3430 


0.24 


0.20 


BS 8752 


G4 


-0.5 


4970 


0.00 


0.00 


BS 6713 


K0.5 lib 


-2.2 


4530 


0.10 


0.08 


BS 6498 


K2 II 


-2.1 


4340 


0.16 


0.14 


BS 2615 


K3 lb 


-1.8 


4080 


0.16 


0.13 


BS 6693 


Ml lb 


-1.6 


3550 


0.26 


0.21 



This paper has been produced using the Royal Astronomical 
Society/Blackwell Science DTjtX style file. 
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Figure Al. Observed variation of the CO spectroscopic index CO^, JW (measured over the original wavelength range proposed by DJW) 
with spectral type and luminosity class in nearby field stars. Compare this figure with the equivalent plot in Appendix A of DJW, which 
employed the spectral atlas of Kleinmann & Hall (1986). The plotted curves are the relations between CO^, JW and T e fj determined by 
DJW for each luminosity class. 
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Figure A2. As for Figure Al, but using the R ~ 3000 stellar spectral atlas of Wallace & Hinkle (1997). 
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Figure A3. Comparison of the CO spectroscopic indices in nearby field stars and in the M100 knot spectra. CO^, JW are the direct 
integrated measurements between rest wavelengths of 2.31—2.40 /im, while CO^ p DW are measurements over the rest wavelength range 
2.293-2.320 (im which are then converted to the equivalent full-range value. The dashed line indicates what would be a 1:1 correspondence. 
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